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Abstract 

We present a microscopic approach for the couphng of cortical activity, as resulting from proper 
dipole currents of pyramidal neurons, to the electromagnetic field in extracellular fluid in presence 
of diffusion and Ohmic conduction. As a result, neural activity becomes represented by a continuous 
neural field equation, while an observation model for electric field potentials is obtained from the 
interaction of cortical dipole currents with charge density in non-resistive extracellular space as 
described by the Nernst-Planck equation. 
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One of the most important problems in the biophysics of neural systems is understand- 
ing the coupling of complex neural network dynamics to the electromagnetic field, that is 
macroscopically measurable (as a feedforward effect) as neural mass potentials, such as local 
field potential (LFP) or electroencephalogram (EEG). Feedback effects, on the other hand, 
have recently been demonstrated via experiments [1]. Thus a theoretical framework for de- 
scribing that coupling — outlined in this RC — is mandatory in clinical, computational and 
cognitive neurosciences, e.g. for treatment of epilepsy {2!, yj or modeling cognition-related 
brain potentials ^|. 

The generators of neural mass potentials are cortical pyramidal neurons (sketched in 
Fig. H]). They exhibit a long dendritic trunk separating mainly excitatory synapses at the 
apical dendritic tree from mainly inhibitory synapses at the basal dendritic tree. When both 
kinds of synapses are simultaneously active, inhibitory synapses generate current sources and 
excitatory synapses current sinks in extracellular space, causing the pyramidal cell to behave 
as a microscopic dipole surrounded by its characteristic electrical field, the dendritic field 
potential (DFP). The densely packed and parallel aligned pyramidal cells form a dipole layer 
whose superimposed currents give rise to LFP and EEG 3, ^• 

Neural mass potentials are most realistically simulated by means of compartmental mod- 
els However, because compartmental models are computationally extremely expansive, 
arge-scale neural networkspreferentially employ point models, based either on conductance 

or population models {q], where mass potentials are estimated either through sums of 
postsynaptic potentials or postsynaptic currents. 

Yet another difficulty is the coupling of the activity of discrete neural networks to the 
continuous electromagnetic field since neural network topology is not embedded into physical 
space as an underlying metric manifold. This problem could be avoided by continuous neural 
networks investigated in neural field theory (NET) Q-Q, and in fact, Q gave first 
accounts for such couplings in NET population models. 

n 

However, both approaches [13] are not concerned with microscopic pyramidal dipole 
currents. Furthermore, extracellular space was assumed to be purely resistive. On the other 
hand, recent research has revealed that diffusion currents, represented by their corresponding 
Warburg impedances [l3|, cannot be neglected in extracellular space as they substantially 



contribute to the characteristic power spectra of neural mass potentials 
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In this RC, we propose a theoretical framework for the microscopic coupling of continuous 
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FIG. 1: Sketch of a cortical pyramidal neuron with extracellular current dipole between spatially 
separated excitatory (open bullet) and inhibitory synapses (filled bullet). Neural in- and outputs 
are indicated by the jagged arrows. The x-axis points toward the scull. Current density j is given 
by dendritic current through cross section area A. 

neural networks (i.e. neural fields) to the electromagnetic field, properly described by dipole 
currents of cortical pyramidal neurons and diffusion effects in extracellular space. As a 
starting point we use a reduced compartment model for a single pyramidal cell and derive 
the evolution law for the activity of a neural network. Additionally, we obtain an expression 
for the dipole current as an observation model. Performing the continuum limit for the 
network yields a neural field equation coupled to the Maxwell equations in extracellular 
fluid. 

We consider populations of neurons, arranged in layers F, {i = 1, . . . , A^). Neurons 
in layers 1 to M should be excitatory, neurons in layers M + 1 to A^ should be inhibitory 



and layer one contains the cortical pyramidal cel. 
arrangement is described by the Amari equation 



s. For a continuous neural network, this 
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TidtUi{x,t) + Ui{x,t) = ^ dx' dt' Wikix,x')si{t - t')fe{uk{x',t')) + hi{x,t) , (1) 

fc=l -^^k J -CO 

where Tj is the characteristic time constant of population i, Ui{x, t) is the neural field activity 
in layer i at time t, Wik{x,x') is the synaptic weight kernel between sites 
Siit) is the postsynaptic impulse response function for layer z, and hi{x,t) is external input 
delivered to the neuron at x in layer Fj. Wave-to-spike conversion is described by the 
sigmoidal activation function fg{u) = + e^*^"^^-*] with activation threshold ^ > 0. 

Herein, we derive an analogous equation for the neural field coupled to the electromagnetic 
field in extracellular space, where the input to layer one turns out as diffusion current hi = 
—kAD-eOxP; with n as some coupling constant, A the current cross section, Einstein's 



diffusion constant [16|], and p the extracellular charge density. This neural field equation 
will be complemented by an observation model for the extracellular dipole current and the 
resulting DFP. 

We describe the zth cortica l py ramidal neuron [Fig. [1] via an electronic equivalent three- 
compartmental model Fig. [2] [171], which is parsimonious to derive a scalar field: one com- 
partment for the apical dendritic tree, another one for the basal dendritic tree, and the third 
for the axon hillock where membrane potential is converted into spike trains. 

Excitatory synapses are represented by the left-most branch, where excitatory postsy- 
naptic potentials (EPSP) at a synapse between a neuron j from layers 1 to M and neuron 
i act as electromotoric forces Efj. These potentials drive excitatory postsynaptic currents 
(EPSC) ifj, essentially consisting of sodium ions, through the cell plasma with resistance 
Rfj from the synapse towards the axon hillock. 

The middle branch describes the inhibitory synapses between a neuron k from layers 
M + 1 to and neuron i. Here, inhibitory postsynaptic potentials (IPSP) Ejf^ provide a 
shortcut between the excitatory branch and the trigger zone, where inhibitory postsynaptic 
currents (IPSC) Iji^ (essentially chloride ions) close the loop between the apical and basal 
dendritic trees. The resistivity of the current paths along the cell plasma is given by Rlj^. 

The cell membrane at the axon hillock itself is represented by the branch at the right 
hand side. Here, a capacitor Ci refiects the temporary storage capacity of the membrane. 
The serial circuit consisting of a battery and a resistor i?^ denotes the Nernst resting 

n 

potential and the leakage conductance of the membrane, respectively [18| . Finally, a gener- 
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FIG. 2: Electronic equivalent circuit for a pyramidal neuron. 

ator of Hodgkin-Huxley spikes [8] (indicated by a "black box" ) is regarded of having infinite 
input impedance. Both, EPSP and IPSP resuh from the interaction of postsynaptic receptor 
kinetics with membrane capacitance of compartments one and two, respectively |l9|. Hence 
these capacitances, omitted in Fig. [21 are already taken into account by Efj,Ejf^ (to what 
we refer as "reduced compartmental model" here). 

The three compartments are coupled through longitudinal resistors, R^,Rf,Rf,Rf 
where , Rf denote the resistivity of the cell plasma j20| and Rf that of extracellular 
space. Yet, in extracellular space not only Ohmic but also diffusion currents are present 
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15|. These are considered by the current source Jf^ connected in parallel to Rf. On 
the other hand, diffusion currents along the somatic resistor Rf are disregarded in adiabatic 
approximation. 

Finally, the membrane voltage at the axon hillock Ui as the dynamical variable and the 
DFP Vi are indicated. The latter drops along the extracellular resistor Rf. For the aim of 
calculation, the mesh currents If (the dendritic current). If (the basal current) and If^ 
(the Hodgkin-Huxley current) are indicated. 
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The circuit in Fig. [2] obeys the following equations: 
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If = (3) 

/f ^ = /f (4) 

/™ = cA + m - E^)/R^ (5) 

Ef, = Rpl + Rfl? + (i?f + )/™ + f/. + (6) 

+ /?P(/°-^P),l<J<P 
El = Rl,ll + iRf + R?)I^'' + U,,l<k<q. (7) 

Here, p is the number of excitatory and q is the number of inhibitory synapses connected to 
neuron i. 

The circuit described by Eqs. ([2] - [7]) shows that the neuron i is likely to fire when 
the excitatory synapses are activated. Then, the Ho dgkin- Huxley current If^^ equals the 
dendritic current Jf. If, by contrast, also the inhibitory synapses are active, the dendritic 
current if follows the shortcut between the apical and basal dendritic trees and only a 
portion could evoke spikes at the trigger zone [Eq. (^]. On the other hand, the large dendritic 
current , diminished by some diffusion current J? , flowing through the extracellular space 
of resistance Rf, gives rise to a large DFP Vi. 

In order to simplify the following derivations, we gauge the resting potential to = 0, 
yielding 

J™ = CA + UJR"^ . (8) 

From Eqs. ([2]), ([6]) and ([H]) we algebraically derive an equation for the extracellular 
current 

p 

^? = Y1 "^^-^ - ^^^^ - ^^^^ + ^^-^^ ' (9) 

with the following electrotonic parameters aij = {Rfj[l + gf{Rf^ + Rf)]}^^ , f3i = Cigf{Rf + 

R?)/[^ + gfiRt + R?)l 7. = 9f{R^ + Rf + R?)/{R''[i + gfiRt + R?)]}, 5. = 9fR?/[i + 

9f{Rt + R?)l^nd gf = EUy^r 

Correspondingly, Eqs. ([3]), (jll), ([7] - [9]) lead to the dynamical law for the membrane 
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potential at axon hillock 

p q 

tA + U, = J2 4 - E ^'^k + ^^^J? , (10) 

3=1 k=l 

where we have introduced the following parameters: time constants: Ti = rjCj[l — g}{Rf + 
Rf) + f3i], excitatory synaptic weights: wfj = r^aij, inhibitory synaptic weights: wjf, = ri/Rli^, 
dendritic diffusion resistance: Ki = ri6i, with fj = R^/[l — gl{Rf + Rf + R^) + 7j-R^] and 
9i = Ylk=i ^l^lk- 

By means of (fTOj) we eliminate the temporal derivative in ([9]), yielding 

p q 

ir = Yl + E ^''^ + ^^^^ + ' (11) 
j=i k=i 

with parameters wfj = aij - f3i/Ti wfj, wj^ = f3i/Ti w}^, = f3i/Ti - 7^, and T]i = 6i - (3i/Ti Ki. 

The change in sign of the inhibitory contribution from Eq. ffTOj) to Eq. ffTTl) has an 
obvious physical interpretation: In (ITO!) . the change of membrane potential Ui and therefore 
the spike rate is enhanced by EPSPs but diminished by IPSPs. On the other hand, the 
dendritic shortcut current If in fllip is large for both, large EPSPs and large IPSPs. 

From Eq. ( ITOi) we eventually get the neural network's dynamics by taking into account 
that postsynaptic potentials are obtained from presynaptic spike rates through temporal 
convolution with postsynaptic impulse response functions, i.e. 

Ef(t)= f sf\t-f)R,{t')dt' (12) 
J —00 

where sf ^(t) are excitatory and inhibitory synaptic impulse response functions, respectively, 
and Rj is the spike rate 

RAt) = h{u,{t))- (13) 

Inserting f|T2|) and f|T3|) into ( ITOj) and performing a continuum limit Ui{t) Ui{x, t) where x G 
Pj replaces the neuron index i, that thereby becomes a population index, entails the Amari 
equation ([1]), with hi = nAj^ and j° = J^/A as diffusion current density through cross 
section A [Fig. [T]. The synaptic weight kernels Wik{x,x') are obtained from the excitatory 
and inhibitory synaptic weights by absorbing their signs after rearrangement to the Amari 
layers. 

In a linear, isotropic medium, currents and gradients are all aligned parallel to the main 
dendritic trunk, indicated by the x-axis in Fig. [1] Therefore, a one-dimensional analysis of 
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neural electrodynamics is sufficient. As further simplification we consider currents of only 
one ionic species (sodium) here, by treating IPSC's as "sodium hole" currents in analogy 
to semiconductor physics. The total current through extracellular fluid is given by the 



Nernst-Planck equation 
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18| 



-D^d^p + aE , (14) 



with Einstein's diffusion constant [16(], charge density p, conductivity a and electric field 
E; the diffusion current in ([H]), j° = —D-^dxP, replaces the input to the ffist layer of the 
Amari equation ([T]): 

TidtUi{x, t) + Ui{x, t) = dx' dt' Wik{x, x)si{t - t')fo{uk{x', t')) - 5i^iKAD^d^p , 

fc=l "^Tfe J -co 

(15) 

where 5i^i = 1(0) if i = l(z 7^ 1) is the Kronecker symbol. Correspondingly, the continuum 
limit of Eq. f lTT]) becomes 



j{x,t) = ^ dx' dt' wik{x,x')si{t -t')fe{uk{x\t')) + ^ui{x,t) - 7]AD^d^p . (16) 
k=i -^^fe J -00 

Moreover, the electric field E is given by the gradient of the DFP V , E = —d^rV and 
the conductivity of the extracellular electrolyte relates to its mobility phj a = pp [18|. In 
addition to the Nernst-Planck equation f[T^ . we have the ffist Maxwell equation dx{eE) = p, 
with permittivity e, and the continuity equation d^j + dtp = reflecting the conservation of 
charge as a result of the Maxwell equations. Computing the divergence of the Nernst-Planck 
equation (|T^ by taking the continuity equation into account, yields, after eliminating E by 
means of f|T^ . dxE by means of the ffist Maxwell equation and a, 

dtp = D^dlp - j9.(lnp) - De^^ - ■ (17) 

p e 

Finally, the DFP is obtained as a solution of ([1] 



dxV = -—(j+Dj,dxp) . (18) 
pp\ ) 

Altogether, Eqs. (fT5|) - (ITS!) describe the coupling of a layered neural field to the elec- 
tromagnetic field in extracellular space, where the pyramidal layer described by the Amari 
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equation ( !T5|) is driven by diffusion currents. Moreover, Eqs. ( !T6|) - ( ITS!) constitute an obser- 
vation model for currents, charge distribution and DFP in extracellular space. The system 
(IT3]) - (ITS]) of nonlinear coupled partial integro-differential equations has to be solved in a 
self-consistent manner, e.g. by means of mean-field techniques. We leave this analysis for 
subsequent research. 

Our model of neural fields coupled to the electromagnetic field can be straightforwardly 
generalized in several directions: (1) Introducing anisotropic media described by conductivity 
and permittivity tensors, respectively, requires full-fledged 3-dimensional calculus. (2) LFP 
and EEG could be obtained as mean-fields from appropriate spatial coarse grainings. (3) 



Neural field effects such as ephaptic interactions 



m 



m 



21| could be phenomenologically 



modeled by volt age- modulation of activation thresholds 6 in the activation function of (ITHj) 



17| 



Since ( !T7|) is a generalized diffusion equation, the present approach is consistent with 



15l |. We expect that our 



the Warburg impedance approach to neural mass potentials |14l . 
theory for the coupling of continuous neural networks to the electromagnetic field in non- 
resistive nervous tissue is significantly important for applications in the clinical and cognitive 



neurosciences, such as research on epilepsy 



3| or cognition-related brain potentials 
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